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Abstract 

The gauge invariant two-point correlator for the gluon field strength tensor 
is analysed by means of the QCD sum rule method. To this end, we make 
use of a relation of this correlator to a two-point function for a quark-gluon 
hybrid in the limit of the quark mass going to infinity. From the sum rules 
a relation between the gluon correlation length and the gluon condensate 
is obtained. We briefly compare our results to recent determinations of the 
field strength correlator on the lattice. 
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1 Introduction 



The gauge invariant non-local gluon field strength correlator plays an important 
role in non-perturbative approaches to QCD fT)-^. It is the basic ingredient in 
the model of the stochastic vacuum (MSV) ||[7| and in the description of high 
energy hadron-hadron scattering f8HTTf . In the spectrum of heavy quark bound 
states it governs the effect of the gluon condensate on the level splittings fl2Hl5f 
and it is useful for the determination of the spin dependent parts in the heavy 
quark potential |T6l,|T7|. 

Its next-to-leading order correction in perturbative QCD has been calculated 



recently by two of the authors |18| . The correlator has also been measured on the 
lattice in pure gauge theory and full QCD using the cooling method [|H||2(| and 
by making the assumptions of the MSV from lattice calculations of the heavy 



quark potentials [21]]. The lattice analyses found that for distances z of the gluon 
field strength larger than roughly 0.4 fm an exponential decaying term dominates 
yielding a correlation length of approximately 0.2 fm. On the other hand the 
short distances are dominated by the perturbative 1/z 4 behaviour. Recently, 
the field strength correlator has also been calculated in the framework of exact 
renormalisation group equations [22]. 

The gauge invariant gluon field strength correlator can be related to a corre- 
lator of a colour singlet current composed of a (fictitious) infinitely heavy octet 
quark and the gluon field strength tensor. This fact has already been employed in 
ref. |18| in order to apply the machinery developed in the Heavy Quark Effective 
Theory]] (HQET) for calculating the perturbative corrections. In this paper we 
again use this relation in order to apply QCD sum rule techniques [|[ to the cor- 
relator in question. The sum rule analysis can be used to estimate the correlation 
length of the field strength correlator using as ingredients the value of the gluon 
condensate and the results for the perturbative calculation. 

Our paper is organised as follows. In the next section we discuss again the 
relation of the field strength correlator and the corresponding heavy quark current 
correlator. In section 3 we set up the different contributions needed for the sum 
rule analysis and in section 4 we present our results together with a comparison 
with recent lattice determinations of the field strength correlator. Finally, in 
section 5, we end with some conclusions and an outlook. 



1 For a review on HQET as well as original references the reader is referred to |2J 
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2 The field strength correlator 



The gauge invariant two-point correlation function of the QCD field strength 
tensor F£ u (x) in the adjoint representation can be defined as 

V, vpa {z) = (OlT^i^CyJPe'^/o^^F^^llO), (2.1) 

where the field strength = d IJ ,A" — d l ,A^ + gf abc A b Al, z = y — x and V denotes 
path ordering of the exponential. In general, the gauge invariant field strength 
correlator could be defined with an arbitrary gauge string connecting the end 
points x and y, but in this work we shall restrict ourselves to a straight line. 
Only for that case the relation to HQET is possible. From the Lorentz structure 
of the field strength correlator it follows that the correlator can be parametrised 
in terms of two scalar functions T>(z 2 ) and T>i(z 2 ) [0: 

£Wx(» = [fiWfiW - 9pa9v P \(^ > {z 2 ) +£>!(> 2 )) 

r _ l dV^z 2 ) 

QppZyZcj QpcrZf/Zp QvpZpZa + QvcjZpZp Qz 2 ' ( ' 

The invariant function T>(z 2 ) can only occur in a non-abelian gauge theory or an 
abelian one with monopoles. In the MSV it is responsible for confinement and 
the formation of a string. 

The correlator T> tlvpa (z) can be related to the correlator of a local, gauge in- 
variant current composed of an infinitely heavy quark field in the octet represen- 
tation, h a (x), and the gluon field strength tensor The current in question 
takes the form {g s h a F^ u ){x). Analogously to HQET the heavy octet-quark field 
is constructed from the field Q a with a finite mass iriq in the limit 

h\x) = lim \ (1+0 e^ x Q a {x) , (2.3) 

with v being the four-velocity of the heavy quark. The propagator of the free 
heavy quark field h^x) in coordinate space is given by 

S(z) = (0\T{h a o (y)h b o (x)}\0) = 8 ab ±6(z°)5(z-^v), (2.4) 

where v° is the zero-component of the velocity. The correlator of the full field can 
be obtained by integrating out only the heavy quark and leaving the expectation 
value with respect to the gauge field: 

(0\T{h a (y)h b (x)}\0) = S(z)(0\Ve 9fabCzT lo dtA ^ +t ^\0}. (2.5) 
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The gauge string is left after the elimination of the heavy quarks from the inter- 
action term of adjoint quarks with the colour potential 



int 



igsfa bc v fl h a (x)A c (x)h b (x) . 



(2.6) 



The physical picture of this result is a heavy quark moving from point x to y with 
a four-velocity v, acquiring a phase proportional to the path-ordered exponential. 
The limit of rriQ — > oo is necessary in order to constrain the heavy quark on a 
straight line and in order to decouple the spin interactions. The same relation 
also holds for quarks in the fundamental representation with the appropriate 
replacements in the exponential. 

The equation ( |2.5| ) allows to establish a relation between the field strength 
correlator (|2.1| ) and the correlator for the colourless heavy quark current. By 
integrating out the heavy degrees of freedom and using (|2.5| ) we arrive at 



We may view the composite operator (g s h a F^ u )(x) as an interpolating field of 
colourless quark gluon hybrids and evaluate V pupcT (z) by introducing these as 
intermediate states in the absorption part of T> fiypa (z). The lowest lying state 
will govern the long-range behaviour and hence the inverse of its energy is the 
correlation length. 

Our next aim is to evaluate this correlator in the framework of QCD sum 
rules |J and in that way obtain information on the correlation length of the 
gluon field strength correlator. For the sum rule analysis it is preferable to work 
with the correlator in momentum space. Thus we define 



where the residual heavy quark momentum is w = vq. Similar to the Lorentz 
decomposition of the coordinate space correlator T> pvp(T {z) into scalar functions 
V(z 2 ) and V 1 (z 2 ), eq. ( |2.2|) , we can write the momentum space correlator as 
follows: 




(0\T{g 2 F^(y)h a (y)F h pa (x)h\x)}\0) 
S(z) D pupa (z) . 



(2.7) 



V,u P a{w) = i J dz e iqz (0\T{g 2 F^ u (y)h a (y)F^ a (x)h b (x)}\0) 



(2.8) 




+ g^pVyVa - g^VuVp - g V pV p V a + gwV^Vp V*(w) . (2.9) 
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The functions T>{w) and T>\(w) are the Fourier transforms of S(z)T>(z 2 ) 
and S(z)T>i(z 2 ) respectively, the function T>*(w) is the Fourier transform of 
S(z)z 2 dV 1 (z 2 )/dz 2 . 

For our purpose of isolating intermediate states of the correlator (|2.8|) a de- 
composition according to an O3 classification is more appropriate than the de- 
composition of eq. (|2.9|) . Since the spin of the heavy quark decouples, we only 
have to consider the gluon spin. The six-component field can be decomposed into 
tensor structures depending on the only two external vectors in the game; the 
four- velocity v p and the polarisation vector of the gluon e p . This leads to the two 
Lorentz structures for the hadronic matrix elements 

(0\(g s F; u h a )(0)\H-) = f( ¥ ,- V ,), (2.10) 
(0\(g s F; u h a )(0)\H + ) = / + £^V, (2.11) 

where H T are hadronic states with the same quantum numbers as the composite 
current. In the rest frame v = 0, the first structure transforms as a 3- vector and 
thus H~ corresponds to a 1~ state whereas the second structure transforms as 
an axialvector and H + corresponds to a 1 + state. 

Through appropriate projections the two quantum numbers can be singled 
out from the correlator V^p^iw). Hence, we define 

V'{w) = g» p v»v a V^ pa (w) = 3{v(w)+V 1 (w)+V,(w)), (2.12) 
V+( w ) = (g^g™ -2g^v°)Vp Vp(7 (w) = § (v{w) +V 1 {w)) . (2.13) 

The Fourier transforms of the functions T>~(w) and T> + {w) are up to the factor 
S(z) the invariant functions T>\\(z 2 ) and D±(z 2 ) respectively which have been 
used in the lattice calculations of refs. |T9|,pO|. 

Under the assumption of quark-hadron duality which is usually made for sum 
rule analyses PJ, we model the correlators by a contribution from the lowest 
lying resonances plus the perturbative continuum above a threshold Sq. Inserting 
the matrix elements and performing the heavy quark phase space integrals one 
obtains 

Z*(«0 = *T' 2 + JdX ; PT(A) , (2.14) 

s o 

where k~ = 1, k + = —2 and E represents the energy of the glue around the heavy 
quark. The spectral densities are defined by p T (A) = I/ttihiI^A + ie) and are 



4 



known at the next-to-leading order [18[]. Explicit expressions will be given in the 
next section. 

After Fourier transformation to coordinate space the above representation 
reads: 

oo 

= \-K\f\ 2 e-^ + J d\p{\)e- a \ z ^S{z), (2.15) 

SO 

where the factorisation of the heavy quark propagator can be seen explicitely. 
The inverse correlation length is found to be given by E. 



3 The sum rules 



The phenomenological side of the sum rules has already been given by eq. ( |2.14j ). 



In this section, we shall present the theoretical side of the sum rules which arises 
from calculating the correlator of eq. (|2.8p in the framework of the operator 
product expansion p|,[25|j. 

In coordinate space the purely perturbative contribution up to the next-to- 
leading order in the strong coupling constant has been calculated in ref. [ |l8|] . 
Here we give the corresponding results in momentum space for V T (w): 

V^iw) = (-w) 3 a [p* + p*L + a (p* + p\L + p^L 2 ) ] , (3.1) 

where a = a s /ir, L = \n(—2w/fi) and the coefficients pfj are given explicitly in 
the appendix. From this result one can immediately calculate the corresponding 
spectral functions: 



p=F(A) = A 3 a 



2A\ 

p^ + a[pT 1 + 2pT 2 \n—j 



(3.2) 



where A has to be greater zero. Essential for the sum rule analysis are the 
contributions coming from the condensates. The correlation function is expanded 
in powers of 1/w corresponding to higher and higher dimensional condensates. 
In our case the dimension three condensate (hh) vanishes since the quark mass 
is infinite. The lowest nonvanishing term is the gluon condensate of dimension 
four: 

7T 2 2ir 2 

Vp F (w) = (aFF) , V+ F (w) = (aFF) . (3.3) 

w w 
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The next condensate contribution would be of dimension six, but we shall neglect 
all higher condensate contributions in this work and restrict ourselves to the gluon 
condensate. 

In order to suppress contributions in the dispersion integral coming from 
higher exited states and from higher dimensional condensates, it is convenient 
to apply a Borel transformation with T being the Borel variable Q. Some 
useful formulae for the Borel transformation are also collected in the appendix. 
For the phenomenological side of the sum rules, eq. fl2.14|) , we then find 



^\n 2 e- E ^' T 



+ / d\p T {\)e 



-X/T 



(3.4) 



s = 



For the perturbative contribution it is convenient to apply the following identity: 



B T V(w) = T A B T (^-) V(w 



from which we obtain 



dw J 



1 2T 
Pfi + a [ v\ + 7T r'(4) + 2 V l 2 In — 
3 fi 



(3.5) 



(3.6) 



where je is Eulers constant and T'(4) = 11 — 67^. The Borel transformed 
expression for the gluon condensate contribution is found to be: 



V 



FF 



7T 



(aFF) 



V+r = 2vr 2 (aFF) 



(3.7) 



After Borel transformation, the correlators satisfy homogeneous renormalisation 
group equations. Thus we can improve the perturbative expressions by resum- 
ming the logarithmic contributions. The perturbative contribution is then ex- 
pressed in terms of the running coupling a(2T): 



1 , a(2T) 



-7f//?i 



a(fi) 



a(2T) 



(3i 



where (3\ = 11/2 — n//3 is the first coefficient of the QCD /3-function. Reexpand- 
ing and comparing with eq. ( |3.6|) , the anomalous dimensions 7^ are found to be 
1\ = ZPm/Pn + Pi> or explicitly 



7i 







3. 



(3.9) 
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Let us note that the correlator T>~(T) which corresponds to the vector interme- 
diate state does not depend on the renormalisation scale /x at this order. 

For the continuum contribution we first evaluate the integral with the general 



formula [p6] which makes the numerical analysis easier: 



r 9X 

/ dX A a_1 ln n — e ~ A / T = T a V 
I A* fro 



rr 



ln fe — 



( so 
da n ~ k V*' T 



(3.10) 



some formulae for the incomplete T-function r(a, x) are given in the appendix. 
We then obtain 



oo , 

X T (T,s ) = JdXpT(X)e- x / T = T'a { pfi r (4, ^ 



+ a 



^ + 2^ 2 in^r(4,|) + 2^ 2 r (4, 1 



and after renormalisation group improvement 



X T (T,s ) 



rp4 



+ a 



a(2T) ^r(4,|) 



/• 1 M ; ^)-2p=r' (4,| 



In the limit So — >■ 0, eq. (|3.12|) agrees with eq. ( p.8|) as it should. 



(3.11) 



(3.12) 



4 Numerical analysis 

After equating the phenomenological and the theoretical part we end up with the 
sum rule 

K T (T) ee - K T in 2 e -EVT = vt f + VMT)-X t (T,s ), (4.1) 

where k~ = 1 and k + = —2. In order to estimate the binding energy we derive 
as an immediate consequence of ( |4.1| ): 







d(l/T) 



In K T 



s_(pj t(t) _ xT(T) 



so, 



Vf F + VT T (T)- X T(T, s ] 



(4.2) 
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The derivative can also be given analytically if we first derive with respect to T 
and then perform the resummation of the logarithms. We thus find 





~ T WJ a(2T) ^H r(5) " r l 5 'T 



+ a 



i>l (i^)-i- (5,|))+2^ (r'(s)-r' (5, | 



(4.3) 



We note that the different signs of the perturbative and non-perturbative terms 
in the 1~ state lead to a stabilisation for the energy sum rule, whereas the equal 
sign in the 1 + state destabilises the sum rule. 

Let us begin our numerical analysis with the case for three light quark flavours. 
As our input parameters we use (aFF) = 0.024±0.012 GeV 4 and A 3fl = 325 MeV. 
In principle, the coupling constant at next-to-leading order could be evaluated 
at any scale \i. As our central value in the numerical analysis we have chosen 
H = 2 GeV. For the energy E~ of the 1~ state we obtain the best stability for 
a continuum threshold Sq = 1.7 GeV in the range T > 0.7 GeV with an energy 
E~ »s 1.4 GeV. To estimate the errors we have varied the scale [i as well as the 
continuum threshold sq. In figure 1 we have displayed the energy E~ as a function 
of the Borel parameter T for // = 1 GeV (dashed lines), 2 GeV (solid lines) and 
4 GeV (dotted lines). The corresponding values of the continuum threshold are 
s = 1.5 ± 0.2 GeV, 1.7 ± 0.2 GeV and 1.9 ± 0.2 GeV respectively. The central 
values have been chosen in order to obtain maximal stability for the sum rule. 

Larger values of so always increase E but at the same time the stability 
region shrinks and goes to smaller values of T. However, even at T = 0.7 GeV 
the influence of the higher resonances expressed through the continuum model 
X~(T) is very large: x~(0-7)/Dp T (0.7) ~ 0.75. For so = 2.1 GeV, we have a small 
stability region around T = 0.65 GeV yielding E~ = 1.6 GeV. Here the influence 
of the continuum model is around tolerable 50%, the perturbative corrections and 
the choice of the renormalisation scale become however more important there. 

Another source of uncertainty is the value of the gluon condensate. For the 
value (aFF) = 0.012 GeV 4 , originally obtained by [0, we find E~ = 1.2 GeV 
at s = 1.5 GeV, whereas for (aFF) = 0.036 GeV 4 we obtain E~ = 1.8 GeV at 
Sq = 2.4 GeV. We therefore conclude from the sum rules for the above mentioned 
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Figure 1: The energy E~ as a function of the Borel-parameter T for three different 
renormalisation scales \i and continuum thresholds Sq. Dashed curves fi = 1 GeV: 
lowest s = 1.3 GeV, middle s = 1.5 GeV, upper s = 1.7 GeV. Solid curves 
/i = 2 GeV: lowest so = 1.5 GeV, middle so = 1.7 GeV, upper so = 1.9 GeV. 
Dotted curves \i = 4 GeV: lowest s = 1.7 GeV, middle s = 1.9 GeV, upper 
s = 2.1 GeV. 



parameters an energy E and a correlation length a of: 

Ez fl = 1.5 ±0.4 GeV and a^ fl = 0.13^o2 fm - ( 4 - 4 ) 

The main sources of uncertainty are the value of the gluon condensate and 
the continuum contribution. Though the perturbative two-loop contributions to 
the sum rule are very large, their influence on the value of E~ is not so dramatic. 



The corrections tend to cancel in the ratio of eq. (|4.2|) . If one determines the 
energy from the sum rule just containing the lowest order perturbation theory 
and chooses as the scale for a s the approximate value of the energy one finds for 
(aFF) = 0.024 GeV 4 the value E~ = 1.9 GeV. 

In a world without light quarks, i.e. rif = 0, the main influence on the 
sum rule is the expected change of the gluon condensate which might increase 
by a factor two to three ||27[| . If we perform an analysis as above, we get for 
A /z = 250 MeV [[§], (aFF) = 0.048 ± 0.024 GeV 4 and s = 2.3 GeV an energy 
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and correlation length of 

E ofl = 1.9±0.5GeV and a^ fl = 0.11+°:o2 fm - ( 4 - 5 ) 

For E + , the energy of the axial vector state, we obtain no stable sum rule. 
Although the expressions for E~ and E + are equal in lowest order perturbation 
theory, higher order perturbative contributions and the gluon condensate lead to 
a splitting in such a way that for the same values of Sq and T the resulting value 
for E~ is higher than that for E + . 



5 Summary and conclusions 

The analysis of the gauge invariant gluon field strength correlator by QCD sum 
rule methods allows to establish a relation between the gluon condensate and 
the correlation length. In order to apply the sum rule technique which consists 
in the comparison of a phenomenological Ansatz with a theoretical expression 
obtained from the operator product expansion we interpret the gluon correlator 
as the correlator of two colour neutral hybrid states composed of a (fictitious) 
heavy quark transforming under the adjoint representation and the gluon field. 
The former serves as the source for the gauge string in the correlator. 

In this approach the decomposition in two invariant functions D + and D~ is 
more appropriate than the decomposition of eq. ( |2.2| ), since D~ receives only con- 
tributions from 1~ and D + from 1 + intermediate states (ignoring the decoupled 
spin of the heavy octet quark). Therefore these functions show simple exponen- 
tial behaviour at large distances and not D and Di. The perturbative expressions 
for D + and D~ are nearly degenerate, but the gluon condensate contributes with 
different sign. It stabilises the sum rule for D~ and destabilises for D + . 

The value of the binding energy for the lowest intermediate 1~ state (the 
inverse correlation length of the correlator) with three flavours is determined to 
be E^jn = l/ajj; ~ 1.5 ± 0.4 GeV and with zero flavours to be Eqj 1 = l/cto/z ~ 
1.9 ± 0.5 GeV. The main sources of uncertainty are the choice of the continuum 
threshold s and the value of the gluon condensate. 

Though we find no stable sum rule for the axial vector state we have from 
the difference of the expressions for the 1~ and 1 + state strong evidence for the 
counterintuitive result that the 1 + state is lighter than the vector state. 

The gauge invariant gluon correlator has been calculated on the lattice using 
the cooling technique ]T^, ^tj. There, the analysis has been made by assuming 
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at large distances an exponential behaviour for the invariant functions D and 
Di, which in light of the present investigation seems less justified than the same 
Ansatz for the functions D~ and D + . The results of the lattice calculation are in 
qualitative, but not quantitative agreement with the sum rule results. The lattice 
researchers find correlation lengths for D and Di, a and a 1; which are degenerate 
within the errors. The computations have been done in quenched QCD and with 
four dynamic flavours of staggered fermions at a bare quark mass of d-m q = 0.01 
where d denotes the lattice spacing. They found |20| : 

E~ = E + = - = 0.90 ± 0.14 GeV for flavours and 
a 

ET = E + = - = 0.58 ±0.10 GeV for 4 flavours. (5.1) 
a 

A preliminary analysis of the lattice data based on an exponential behaviour 
for D + and D~ |29|] leaves the values essentially unchanged but indicates a split- 
ting of E + and E~ in the same direction as proposed by the sum rules! The 
reader should also note the increase of the correlation length from zero to four 
flavours which is predicted by the sum rules as well where it is mainly due to the 
decrease of the gluon condensate. 

In another approach |H[ the exponential behaviour of the functions D + and 
D* = z 2 d/dz 2 Di for quenched QCD could be extracted by analysing field in- 
sertions into a Wilson loop and assuming factorisation as in the model of the 
stochastic vacuum || [/J. The resulting values for the correlation lengths are 
smaller than those of the direct lattice calculations [21] and thus compare more 
favourably with our results: 

E + = — = 1.64 GeV and E* = — = 1.04 GeV for flavours. (5.2) 

a+ a* K ' 

The sum rule analysis shows that the state investigated here namely a gluon 
confined by an octet source has a much higher energy than the corresponding 
state in HQET. A similar analysis of a light quark bound by a source in the 
fundamental representation |30| yielded an energy which is by a factor 2 to 4 
smaller. This is to be expected on general grounds p7| since the case treated 
here is nearer to a glueball than to a heavy meson. 
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Appendix 

The theoretical expression for the perturbative correlator up to next-to-leading 
order is given by: 

V* T (w) = (-w) 3 [a(pf + p*L) + a 2 04, + p*L + p* 2 L 2 ) 

where L = ln(— 2w/fj) and the pfj have the values 





40 


Pw = 


~3~ 


Pu = 


- 16 




2839 


P20 = 


9 



692 _ 2 104 
P21 = - — - 167r + ~ n f 

8 

P22 = 44 - -n f 

. 128 
P10 = - — 

Pti = 32 

+ 5684 2 1 _, /oX /848 8tt 2 \ 

Plo = -^--44vr 2 + l92C(3)+^— + — j n f 

+ 1072 2 208 

P21 = ~^~ + 32n ~~ n s 

, 16 
p+ = -A0 + —n f . 

For the convenience of the reader we also give the definition of the Borel 
transformation and some useful formulae: 

Bt — iim £^(ff, r = ^>oaxed 

-TO,n->oo 1 (n + 1) \dw ) n 



g 1 = 1 C -E/T 

T {E-w-ie) a ~ T^T"- 1 
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Below, we have collected some formulae for the incomplete Gamma function 
which are helpful for the numerical analysis of the sum rules: 

T{a,x) = / e'H^dt 

F(n,x) = T(n)e- x J2 ^, n = 1,2,... 

fc=0 

d 

T'(a,x) = — — r(a,x) 

V(a) - T'(a, x) = f e-H a - l \nt dt 

Jo 

r'(4) - T'(4, a;) = 11 - 67^ - 6r(0, x) 

-e~ x (ll + 5x + x 2 + (6 + 6x + 3a; 2 + x 3 ) lnx) 

r'(5) -r'(5,x) = 50 -24 7 £-24r(0,x) 

- e~ x ( 50 + 26a; + 7x 2 + x 3 + (24 + 24a; + 12a; 2 + 4a; 3 + a; 4 ) lnx) . 
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